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The  hybrid  system  of  a  photovoltaic  solar-cell  module  and  a  thermoelectric  module  has  drawn  consider¬ 
able  attention  recently.  Here,  unlike  most  of  the  previously  reported  work,  we  present  an  analysis  that 
intimately  couples  the  nodal  temperatures  of  the  system  and  power  outputs  of  both  modules.  Within 
the  range  of  selected  solar  constants,  we  have  obtained  efficiencies  for  hybrid  systems  in  comparison 
with  those  of  photovoltaic  solar  cells  alone.  Our  findings  indicate  that,  for  parametric  values  chosen 
herein,  the  former  are  generally  lower  than  the  latter,  unless  thermal  electric  materials  of  low-valued 
thermal  conductivities  and  high-valued  Seebeck  coefficients  can  be  found  or  fabricated. 

©  2014  Elsevier  Ltd.  All  rights  reserved. 


1.  Introduction 

Photovoltaic  solar  cells  are  devices  that  convert  solar  energy 
into  electricity. 

During  the  conversion  process,  a  portion  of  the  impinging  radi¬ 
ation  is  absorbed  by  the  cells,  causing  the  rise  of  the  photovoltaic 
temperature.  It  is  known  1  ]  that,  for  photovoltaic  cells,  the  output 
power  decreases  as  the  temperature  increases.  Therefore,  if  we  are 
able  to  design  a  hybrid  system  consisting  of  a  solar-cell  unit  and  a 
thermoelectric  module,  allowing  the  latter  to  dissipate  the  ab¬ 
sorbed  energy  in  the  former,  and  at  the  same  time  to  reduce  the 
photovoltaic  temperature,  there  may  be  some  hope  to  produce  a 
larger  amount  of  the  output  electric  power  than  the  solar-cell 
module  alone. 

During  recent  decades,  advantages  of  thermoelectric  devices 
have  gradually  been  realized  because  these  devices  are  capable 
of  directly  converting  the  thermal  energy  into  the  electricity  within 
small  sizes.  Characteristics  of  thermoelectric  modules  alone  were 
investigated  [2-4].  In  [5],  researchers  started  focusing  on  heat  dis¬ 
sipation  of  a  hybrid  system.  A  couple  of  years  later,  van  Sark  [6] 
adopted  a  new  type  of  thermoelectric  material  when  designing 
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the  device.  A  hybrid  system  with  water  cooling  was  also  studied 
[7].  Using  heat  flux  concentrators,  Vorobiev  et  al.  [8]  developed  a 
hybrid  system  that  yielded  high  efficiencies.  In  [9],  evacuated-tube 
heat-pipe  solar  collectors  were  studied,  with  agreeable  compari¬ 
sons  of  theoretical  results  and  experimental  data.  A  computational 
study  [10  using  the  finite  element  method  was  conducted  for  a 
hybrid  system  with  a  collector  panel  inserted  between  the  photo¬ 
voltaic  cell  and  the  thermoelectric  module,  yielding  high  thermal 
efficiencies.  Finally,  in  [11  ],  the  thermoelectric  module  was  compu¬ 
tationally  investigated,  with  the  solar  flux  input  as  one  of  the 
boundary  conditions. 

Here  in  this  analysis  we  have  placed  the  emphasis  on  the  cou¬ 
pling  between  discretized  nodal  temperatures  in  the  heat  conduc¬ 
tion  phenomenon  and  the  electrical  power  outputs  produced  by 
the  hybrid  system.  In  other  words,  nodal  temperatures  are  func¬ 
tions  of  power  outputs,  and  vice  versa.  Therefore,  T(i  )*  Pte>  and 
Ppv  are  strongly  inter-related.  The  thermoelectric  efficiency,  rjte ,  is 
not  assumed  to  be  the  maximum  value.  These  variables  must  be 
computed  simultaneously.  Such  a  coupling  approach  constitutes 
a  nonlinear  challenge,  and  is  presented  below. 

2.  Modeling  and  formulations 

The  system  schematic  of  an  assembly  consisting  of  a  photovol¬ 
taic  solar  unit  and  a  slab  of  thermoelectric  material  is  shown  in 
Fig.  1.  The  sunlight  impinges  upon  the  top  surface  of  the  assembly. 
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Nomenclature 

4C 

area  of  the  hybrid  system  (m2) 

pj 

energy  of  Joule  heating  (W) 

bi  i 

Biot  number  defined  as  hiAytelkte  for  bottom  face 

Ppv 

power  output  of  photovoltaic  solar  cells  (W) 

bis 

Biot  number  defined  as  h5Ayglkg  for  top  face 

Pte 

power  output  of  thermoelectric  module  (W) 

c 

speed  of  light,  3  x  108(m/s) 

P hybrid 

power  output  of  the  hybrid  system  (W) 

di,d2 

empirical  constants  given  in  Eq.  (6) 

Rr 

radiative  recombination  rate  at  zero  quasi-Fermi  level 

e 

elementary  charge  1.6  x  10“19(C) 

(same  unit  as  Ngap) 

E 

energy  of  a  solar  photon  (eV) 

R,r 

external  and  internal  resistances  of  thermoelectric  mod¬ 

Eg 

bandgap  of  silicon  material  (eV) 

ule  circuit  (Q) 

h 

Planck  constant  6.63  x  10_34(J  s) 

S 

Seebeck  coefficient  of  thermoelectric  module  (pV/K) 

hi 

heat  transfer  coefficient  for  the  bottom  surface  (W/ 

Sc 

solar  irradiance  (solar  constant)  (W/m2) 

m2  K) 

Ti . T5 

nodal  temperatures  at  node  1,. .  .,5  (K) 

h5 

heat  transfer  coefficient  for  the  top  surface  (W/m2  K) 

Pair 

temperature  of  the  ambient  air  (K) 

I 

electrical  current  in  thermoelectric  module  (A) 

Tsurr 

temperature  of  the  surroundings  (K) 

Ipv 

electrical  current  in  the  solar  cell  (A) 

Vte 

voltage  of  the  thermoelectric  module  circuit  (V) 

kte 

thermal  conductivity  of  the  thermoelectric  module  (W/ 

Vpv 

external  voltage  of  the  solar  cell  (V) 

m  K) 

Ayte 

half  thickness  of  the  thermoelectric  module  (m) 

1<sub 

thermal  conductivity  of  the  substrate  in  solar  cells  (W/ 

Ay  sub 

half  thickness  of  the  substrate  (m) 

m  K) 

Ayg 

half  thickness  of  the  glass  layer  (m) 

kg 

thermal  conductivity  of  the  glass  layer  in  solar  cells  (W/ 

m  K) 

Greek  symbols 

n 

number  of  p-n  pairs  in  the  thermoelectric  module 

a,  p 

constants  related  to  bandgap  and  given  in  Eq.  (8) 

Ne 

Ngap  per  unit  electron  volt 

s 

emissivity  of  the  protective  glass  layer 

bigap 

number  of  solar  photons  per  unit  time  per  unit  area, 

P 

reflectivity  of  the  protective  glass  layer 

carrying  energy  larger  than  the  bandgap 

G 

Stefan-Boltzmann  constant  (5.67  x  10-8  W/(m2  I<4)) 

Qs 

heat  flow  rate  entering  the  control  volume  (W) 

Q  hybrid 

efficiency  of  the  hybrid  system 

Qn 

heat  flow  rate  leaving  the  control  volume  (W) 

A  sheet  of  photovoltaic  thin  film  is  imbedded  inside  the  unit  filled 
with  the  protective  glass  and  the  substrate  material,  respectively, 
above  and  below  the  film.  Underneath  the  solar-cell  unit,  a  ther¬ 
moelectric  module  is  attached.  Both  of  the  bottom  surface  and 
the  top  surface  of  the  assembly  are  cooled  by  air  flows. 

The  system  is  assumed  to  be  in  steady  state,  and  the  heat  trans¬ 
fer  phenomenon  takes  place  in  the  one-dimensional  manner,  with 
all  variables  being  functions  of  the  y  coordinate  only.  Based  on  the 
first  law  of  thermodynamics,  we  can  derive  five  nonlinear  algebraic 
equations  governing  Tlt  T2,  I3,  T4,  and  T5,  as  shown  in  Fig.  1.  Even 
though  we  are  aware  that  heat  fluxes  flow  downward  in  real-life 
situations,  we  will  proceed  to  assume  that  they  flow  upward  along 
the  positive  y  direction  to  avoid  sign  confusions. 

For  Tlt 

AMTair  -  Ti)  =Ackte(T !  -  72)/Ayte  +  Pte(l),  (1) 

which  can  be  simplified  and  rearranged,  with  final  coefficients 
listed  in  Table  1.  The  power  output,  Pte(  1),  can  be  expressed  by 
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P*(1)  =  V(1)J-P,(1),  (2a) 

in  which  the  Seebeck  electrical  current,  /,  is  needed.  The  relation¬ 
ship  between  the  current  and  the  Seebeck  voltage  can  be  expressed 
by 

I  =  (V  (1 )  +  V(2)  +  V(3))/(R  +  r(l)  +  r(2)  +  r(3)),  (2b) 

Pj(l)  =  /2r(l),  (2c) 

where  V(l)  =  nS(T,.5  -Ti),  r,.5  =  (T,  +  T2)/2,  (2d) 

and  r(l)  =  e,(r0  +  ^T,).  (2e) 


Values  of  e^,  r0,  and  fJ  are  given  in  Appendix  A.  Definitions  of  volt¬ 
ages,  internal  electrical  resistance,  and  midpoint  nodal  tempera¬ 
tures  at  other  nodes  can  also  be  found  in  Appendix  A.  The 
constant  n  is  the  number  of  p-n  semi-conductor  pairs,  and  S  is 
the  Seebeck  coefficient.  It  can  be  observed  that,  via  Eqs.  (1,  2a-e), 
the  electrical  current  and  the  nodal  temperatures  are  strongly 
coupled. 

For  the  derivation  of  the  governing  equation  for  I2,  the  control 
volume,  over  which  the  energy  conservation  is  taken,  is  shown  in 
the  dashed  line  in  Fig.  2.  The  nodal  temperature  T2  is  situated  at 
the  center  of  the  control  volume.  From  south,  the  heat  flux  enters 
the  control  volume.  It  leaves  the  control  volume  to  north.  The  Joule 
heating  can  be  treated  as  a  heat  sink  inside  the  control  volume,  and 
does  not  explicitly  appear  in  the  energy-conservation  equation.  In 
addition,  there  is  power  output  due  to  Seebeck  thermoelectric  ef¬ 
fects.  These  quantities  are  designated,  respectively,  as  qs,  qn ,  and 
Pte.  Consequently,  according  to  the  flux  directions  indicated  by  ar¬ 
rows,  we  obtain 


Qs  -  Qn  -  Pte(2)  =  0, 


(3) 


Fig.  1.  The  system  schematic  representing  the  hybrid  system  of  photovoltaic  solar 
cell  and  thermoelectric  module. 


where 

Qs  =  Ackte(h  -  T2)/Ay[e,  and  qn  =  Ack[e(T2  -  T3)/Ayte. 
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Table  1 

The  matrix  coefficients  for  the  nonlinear  set  of  5  equations  governing  nodal  temperatures,  power  output  of  the  thermoelectric  module,  and  power  output  of  the  photovoltaic  solar 
cell,  along  with  the  constant  vector. 
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Fig.  2.  A  typical  control  volume  containing  nodal  temperature,  T2.  It  is  noted  that  P,  and  Pte  should  bear  the  same  sign. 


The  thermoelectric  power  output,  Pte( 2),  can  be  obtained  in  a 
similar  manner  as  Eq.  (2a). 

In  the  derivation  of  the  governing  equation  for  T3,  which  is  lo¬ 
cated  at  the  interface  between  the  thermoelectric  material  and 
the  photovoltaic  module,  we  idealize  the  thermal  contact  resis¬ 
tance  to  be  zero.  After  straightforward  algebraic  manipulations, 
we  obtain 

a(3, 2  )T2  +  a(3, 3)T3  +  a(3, 4)T4  =  b(  3),  (4) 

where  the  coefficients  are  listed  in  Table  1. 

The  derivation  of  the  governing  equation  for  T4  is  similar  to  that 
for  I3.  The  final  equation  can  also  be  cast  into  the  standard  form  as 

a(4,3)T3  +  a(4,4)T4  +  a(4, 5  )TS  =  b(  4).  (5) 

It  is  worth  noting  that  the  electrical  power  delivered  by  the  photo¬ 
voltaic  film  can  be  related  to  the  nodal  temperature  by 

Ppv  —  A:(dl  +  d2T4),  (6) 

where  diand  d2  are  constants,  whose  values  depend  on  the  imping¬ 
ing  solar  constants,  and  can  be  computed  a  priori  independently  of 
the  system  characteristics.  Such  a  linear  approximation  greatly  sim¬ 
plifies  the  analysis,  enabling  us  to  officially  couple  the  nodal  tem¬ 
peratures  and  the  power  outputs. 

The  details  are  briefly  given  below. 

First,  the  number  of  solar  photons  per  unit  time  per  unit  area, 
carrying  energy  (that  is  larger  than  the  bandgap)  per  electron  volt 
can  be  given  as 

NE  =  {AM\.5)(hc/E3),  (7) 

where  AMI. 5  stands  for  NREL’s  (National  Renewable  Energy  Labora¬ 
tory)  data  for  the  solar  spectrum  at  Sc=  1000  W/m2,  with  its  unit 
being  W/m2  nm;  h  represents  the  Planck  constant, 
6.63  x  1 0— 34  J  s;  the  symbol,  c,  is  the  speed  of  light,  3  x  108  m/s. 


Next,  according  to  [12],  the  bandgap,  Eg,  decreases  with  increas¬ 
ing  temperatures,  and  can  be  written  as 

Eg(T)=Eg(0)-uT2/(T  +  p),  (8) 

where  Eg{ 0)  =  1.1557  eV,  a  =  7.021  x  1(T4  eV/K,  and  p  =  1108  K. 

Only  photons  that  carry  the  energies  greater  than  the  bandgap 
are  capable  of  exciting  the  electrons  located  in  the  valence  band, 
and  force  them  to  jump  to  the  conduction  band.  The  number  of 
such  photons  can  be  expressed  by 

Ngap  =  NEdE  =  /  (AMI  .5)  dE  =  function  of  T.  (9) 

A  portion  of  the  excited  electrons  may  return  to  the  valence  band, 
and  recombine  with  holes.  For  such  recombination  activities,  the 
rate  can  be  derived  as 

271  r  E2dE 

c2h\lEg  explE/(kBT)}  -  1  ’  liuj 

thus  inviting  Eq.  (8)  to  enter  the  problem  [13].  The  electrical  current 
is  generated  from  the  electron-hole  pairs  that  do  not  recombine. 
Note  that  the  only  recombination  mechanism  considered  in  this 
analysis  is  caused  by  radiation.  Under  this  idealization,  we  can  write 

Ipv  —  PNgap  —  6  Rr{6xp[eVpV/ (kpT)]  —  1},  (11) 

where  the  internal  series  resistance  and  the  shunt  resistance  of  so¬ 
lar  cells  have  been  neglected.  Finally,  the  output  power  of  the  solar 
cell  can  be  expressed  by 

Ppv  —  V pvlpv  —  function  of  1/p^.  (1^) 

Under  the  assumption  that  solar  cells  deliver  the  maximum  power 
output,  we  can  numerically  find  the  corresponding  maximum 
power  output  from  Eq.  (12).  The  laborious  procedure  of  finding  dA 
and  d2  can  be  briefly  described  as  (a)  give  a  value  of  T;  (b)  compute 
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Table  2 

Constant  values  of  dx  and  d2  for  various  values  of  Sc. 


Sc(W/m2) 

di 

d2 

800 

440.070 

-0.579 

1000 

510.337 

-0.592 

2000 

1064.852 

-1.287 

3000 

1532.460 

-1.661 

4000 

2043.820 

-2.173 

5000 

2555.176 

-2.677 

6000 

3066.575 

-3.173 

Fig.  3.  Power  output  of  photovoltaic  solar  cell  versus  the  temperature  of  the  thin- 
film  solar  cell,  parametrized  in  solar  constants. 


Eg  according  to  Eq.  (8);  (c)  find  maximum  Ppv;  finally  (d)  use  the 
curve  fitting  to  obtain  and  d2.  Table  2  lists  their  values  for 
Sc  =  800-6000  W/m2. 

In  Fig.  3,  a  family  of  nearly  linear  lines  of  Ppv  versus  the  photo¬ 
voltaic  temperature  are  plotted,  parameterized  in  six  various  val¬ 
ues  of  Sc  for  silicon  solar  cells.  Finally,  it  is  also  noted  that  we 
assume  the  film  to  be  so  thin  like  a  sheet  (of  a  few  nanometers) 
that  its  thickness  is  negligible  in  comparison  with  those  of  thermo¬ 
electric  module,  the  substrate,  and  the  protective  glass  layer. 

For  the  derivation  of  the  governing  equation  for  T5,  we  take  the 
energy  balance  over  the  half-grid-interval  control  volume  between 
x4.5  (the  average  of  x4  and  x5  and  the  top  surface.  This  balance  dic¬ 
tates  that  the  northbound  heat  conduction  due  to  the  temperature 
gradient  between  node  4  and  node  5  should  be  equal  to  the  energy 
leaving  the  sheet  (thus,  energy  components  entering  the  sheet 
should  be  considered  negative),  eventually  yielding 

I4  -  (1  +  bi5)T5  =  c5Aysub/ksub,  (13) 

where  c5  includes  energy  components  of  radiation  and  convection 
entering  and  exiting  above  the  top  surface,  and  is  defined  as 


Cs  =  qr,„et  +  (p  -  l)Sc  -  h5Tair.  (14) 

In  Eq.  (14),  the  first  term  denotes  the  net  radiation  loss  from  the 
system  to  the  surroundings,  and  is  written  as 

qrm=€a(T45-T4sun).  (15) 

Equation  (15)  holds  under  the  assumption  that  the  area  of  the 
absorbing  top  surface  of  the  hybrid  system  is  much  smaller  than 
that  of  the  surroundings,  such  as  the  sky.  In  the  present  analysis, 
this  assumption  is  certainly  valid.  In  addition,  it  should  be  made 
clear  that  the  temperature  of  the  surroundings,  such  as  clouds  in 
the  sky,  may  not  necessarily  be  the  same  as  that  of  the  ambient 
air  near  the  earth  ground. 

Equations  (1),  (2a)-(2e),  (3)-(15)  constitute  a  set  of  nonlinear 
equations  governing  the  characteristics  of  the  assembly  system. 
When  the  solution  converges,  the  efficiency  of  the  hybrid  system 
can  be  defined  and  computed  by 

^1  hybrid  =  ^ hybrid/^  =  [Pte(l)  +  Pfe(2)  +Pte(3)  +Ppv\/SC.  (16) 

The  solution  procedure  is  described  in  the  next  section. 

3.  Computational  procedure 

Equations  (1),  (2a)-(2e),  (3)-(16)  can  be  solved  using  the  stan¬ 
dard  Newton-Raphson  (NR)  method  [14].  However,  we  are  con¬ 
cerned  that,  after  interspersing  linearized  terms  (related  to 
photovoltaic  and  thermoelectric  behaviors)  with  linear  terms  dis¬ 
cretized  from  the  heat  conduction,  we  may  encounter  difficulties 
in  debugging  the  code  in  the  event  of  divergences.  Therefore,  we 
have  opted  the  adoption  of  a  modified  version  of  NR  method,  in 
which  two  stages  of  computations  are  conducted.  In  the  first  stage, 
nonlinear  terms  are  left  in  the  right-hand  side  of  equations.  For  the 
sake  of  convergence  safety,  the  linearized  solution  obtained  by 
using  a  direct  linear  solver  is  preferably  under-relaxed.  In  the  sec¬ 
ond  stage,  all  nonlinear  terms  are  computed  using  the  newly  iter¬ 
ated  and  under-relaxed  values.  The  procedure  then  repeats. 

It  is  understood  that  this  two-staged  under-relaxed  NR  version 
is  clearly  not  as  robust  as  the  standard  NR  version,  and  its  conver¬ 
gence  rate  is  also  slower  than  the  standard  NR  counterpart.  The 
advantage  of  this  modified  version,  however,  is  twofold:  (1 )  a  great 
amount  of  labor  associated  with  algebraic  manipulations  and  cod¬ 
ing  is  saved,  and  (2)  debugging  the  code  is  easier  when  divergence 
happens.  Regardless  of  choices,  the  two  solutions  computed  using 
these  two  methods  should  become  identical  upon  solution  conver¬ 
gence.  The  concept  of  this  mathematical  modeling  is  similar  to  that 
of  the  compartment  modeling  widely  used  in  ecological  sciences 
[15]. 

Fig.  4  shows  the  flow  chart  of  the  computer  code.  In  the  first 
stage,  nodal  temperatures  are  computed  iteratively,  with  power 
outputs  of  the  hybrid  system  either  guessed  or  previously  iterated. 
Iterations  are  needed  because  the  radiation  term,  T4,  exists  in  Eqs. 
(13)-(15).  In  the  second  stage,  these  more  accurate  nodal  temper¬ 
atures  are  used  to  update  values  of  internal  electrical  resistances  in 


Fig.  4.  Flowchart  of  the  algorithm  that  is  used  to  iteratively  compute  nodal  temperatures,  power  output  of  the  thermoelectric  module,  and  power  output  of  the  photovoltaic 
solar  cell. 
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Fig.  5.  Comparison  of  the  temperature  of  the  thin-film  solar  cell  for  two  cases:  (a) 
photovoltaic  solar  cell  only,  and  (b)  the  hybrid  system. 


Fig.  6.  Efficiency  versus  external  electrical  resistance.  The  R  value  that  corresponds 
to  the  maximum  efficiency  is  chosen  for  computations. 


Fig.  7.  Comparison  of  efficiencies  for  two  cases:  (a)  photovoltaic  solar  cell  only 
(rj  =  PpvISc),  and  (b)  the  hybrid  system. 


Fig.  8.  (a)  Efficiency  of  the  hybrid  system  versus  the  thermoelectric  thermal 
conductivity,  and  (b)  power  output  of  the  hybrid  system  versus  the  thermoelectric 
thermal  conductivity.  Both  cases  are  parametrized  in  Sc. 

the  thermoelectric  module  and  these  power  outputs.  For  all  runs, 
the  global  energy  balance  has  been  checked  to  ensure  the  correct¬ 
ness  of  the  code. 

4.  Results  and  discussion 

Fig.  5  shows  the  photovoltaic  temperature,  T4,  versus  the  solar 
constant,  Sc  for  both  cases  of  (a)  the  photovoltaic  thin  film  only 
and  (b)  the  hybrid  assembly  of  photovoltaic  thin  film  and  thermo¬ 
electric  module.  Both  temperatures  are  almost  linear  functions  of 
Sc.  They  are,  however,  associated  with  different  slopes.  It  is  ob¬ 
served  that  values  of  the  hybrid  case  are  overall  higher  than  those 
of  the  thin  film  case,  suggesting  that  the  addition  of  the  thermo¬ 
electric  module  blocks  the  cooling  effect.  As  Sc  increases,  the  differ¬ 
ence  in  two  temperatures  also  increases. 

To  choose  an  appropriate  value  of  external  electrical  resistance, 
R ,  we  have  first  evaluated  the  efficiencies  of  the  hybrid  system  cor¬ 
responding  to  various  values  of  given  R.  Then  the  value  that  yields 
the  maximum  efficiency  is  chosen,  as  shown  in  Fig.  6.  For  example, 
for  a  typical  run  with  the  following  data: 

5c  =  1000  W/m2,  kte  =  0.48  W/(m  K),  n  =  100,  Lte  =  0.003  m, 
LSUb  =  0.003  m,  Lgiass  =  0.003  m,  r  =  3.9973  Q  (computed),  and 
R  =  3.9925  Cl, 

the  efficiency  of  the  hybrid  system  is  equal  to  0.3311. 
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In  Fig.  7,  we  present  the  efficiency,  rj,  versus  the  solar  constant, 
Sc  for  both  cases  of  (a)  the  photovoltaic  thin  film  only  and  (b)  the 
PV-TE  hybrid  assembly.  Both  curves  exhibit  the  trend  that  the  effi¬ 
ciency  increases  as  the  solar  constant  increases  as  expected.  With¬ 
in  the  range  of  the  investigated  Sc  values,  the  efficiency  of  the 
hybrid  system  is  lower  than  that  of  the  photovoltaic  solar  cell 
alone.  This  finding  offers  a  warning  guide  for  the  industrial  manu¬ 
facturers  who  may  be  interested  in  fabricating  such  hybrid  sys¬ 
tems,  with  the  hope  that  the  hybrid  systems  will  yield  higher 
efficiencies. 

Finally,  Fig.  8a  and  b  depict  the  efficiency,  hybrid  and  the  com¬ 
bined  power  output,  Pbybnd  versus  the  thermal  conductivity  of  ther¬ 
moelectric  module,  kte .  We  note  that  both  rjhybrid  and  Phybrid 
decrease  as  k  increases  for  kte  <  0.25.  Beyond  this  kte  value,  both 
the  efficiency  and  the  power  output  remain  fairly  unchanged.  An¬ 
other  noteworthy  feature  is  that  the  sensitivity  of  r\bybrid  on  kte  is 
quite  high  within  the  range  of  low  kte  values.  If  thermoelectric 
materials  having  low  kte  values  can,  indeed,  be  fabricated,  it  ap¬ 
pears  possible  for  the  industry  to  design  a  hybrid  system  that  is 
capable  of  yielding  high  values  of  thermal  efficiencies. 

With  the  advances  on  thermoelectric  research,  such  as  the 
enhancement  of  power  factor  (electrical  conductivity  times  square 
of  Seebeck  coefficient)  through  optimizing  transport  distribution 
functions  [16]  or  nanostructure  designs  [17,18],  and  the  reduction 
of  thermal  conductivity  by  nanostructure  scattering  [17,19]  or 
strain  engineering  [20],  the  efficiency  of  thermoelectric  energy 
conversion  can  be  greatly  increased.  The  PV-TE  hybrid  system 
can  thus  be  potentially  advantageous  due  to  such  significant 
improvements  of  thermoelectric  materials. 

It  is  worth  mentioning  as  a  tender  reminder  that,  in  our  analy¬ 
sis,  the  power  output  of  the  solar-cell  module  is  approximated  by  a 
linear  function  of  the  solar-cell  temperature,  with  constants,  di  and 
d2,  comprehensively  computed  according  to  Eqs.  (7)-(12).  The  ben¬ 
efit  of  this  approximation  is  that  the  power  output  is  reliably  cou¬ 
pled  with  the  nodal  temperatures  of  the  PV-TE  hybrid  system. 

5.  Conclusions 

In  this  study,  we  have  emphasized  on  coupling  the  temperature 
field  with  the  power  outputs  of  the  hybrid  system  that  consists  of  a 
photovoltaic  solar-cell  module  and  a  thermoelectric  module.  The 
efficiency  of  the  thermoelectric  module  is  directly  computed  based 
on  thermoelectric  theories,  and  the  value  of  the  external  electrical 
resistance  is  taken  to  be  the  one  corresponding  to  the  maximum 
value  of  efficiency.  Primary  findings  can  be  summarized  below: 

1.  According  to  our  findings  within  the  parametric  values  chosen, 
the  thermal  efficiency  of  the  hybrid  system  is  generally  lower 
than  that  of  the  solar-cell  module  alone.  If  we  are  able  to  fabri¬ 
cate  thermoelectric  materials  with  low-valued  thermal  conduc¬ 
tivities  and  high-valued  Seebeck  coefficients,  it  may  be  possible 
for  us  to  achieve  rjhybrid  >  rjpv  only. 

2.  The  power  output  of  the  solar  cell  module  can  be  approximated 
a  priori  by  a  linear  function  of  the  temperature  7(4),  which  is 
then  coupled  with  other  nodal  temperatures.  Such  an  approxi¬ 
mation  greatly  facilitates  the  thermal  analysis. 

3.  A  numerical  procedure  has  been  developed  to  couple  the  power 
outputs  of  both  the  photovoltaic  solar-cell  module  and  the  ther¬ 
moelectric  module  with  the  temperature  distribution  of  the 
hybrid  system.  It  serves  as  a  tool  for  the  future  research  in 
designing  PV-TE  hybrid  assemblies. 
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Appendix  A 


e-i 

2.5  xlO4 

e2 

5.0  xlO4 

e3 

2.5  xlO4 

r0 

l.OxlO”5  Q 

P 

1.0  x  10-7  Q/K 

r(l) 

e-i  (r0+/?T,) 

internal  electrical  resistance 

HU 

nS(Tl5-T,) 

voltage 

r-i.5 

(Ti+T2)/2 

midpoint  nodal  temperature 

r(  2) 

e2(ro+/?r,) 

internal  electrical  resistance 

V(2) 

nS(T2.  s-TYs) 

voltage 

T2.5 

(r2+r3)/2 

midpoint  nodal  temperature 

r(  3) 

e3(r0+Ph) 

internal  electrical  resistance 

V(3) 

nS(T3-T2.5) 

voltage 
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